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CD , ABSTRACT 

Q 

We investigate the possibility of an intrinsic Baldwin Effect (i.e., nonlinear 
emission-line response to continuum variations) in the broad H/3 emission line of 
>■ I the active galaxy NGC 5548 using cross-correlation techniques to remove light 

travel-time effects from the data. We find a nonlinear relationship between the 
^ ' H/3 emission line and continuum fluxes that is in good agreement with theoretical 

predictions. We suggest that similar analysis of multiple lines might provide a 
^ ■ useful diagnostic of physical conditions in the broad-line region. 

Ph| Subject headings: galaxies: active — galaxies: individual (NGC 5548) — galaxies: 

O ' nuclei — galaxies: Seyfert 

^' 

1. Introduction 

■ Correlations between the continuum and emission-line properties of active galactic nu- 
clei (AGNs) afford potentially important probes of the structure and physical conditions of 
the broad-line region (BLR) and of the unobservable ionizing continuum in these sources. 
Probably the best-established correlation between continuum and emission-line properties is 
the "Baldwin Effect" (Baldwin 1977), an inverse, non-linear correlation between the equiv- 
alent width of the broad UV emission lines, most notably ClvA1549, and the luminosity of 
the underlying continuum (e.g., Osmer, Porter, & Green 1994). The Baldwin relationship 
can be parameterized as a simple power law of the form 

-^line OC -^cont' (-'-) 
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with index a ^ 1. 

Kinney, Rivolo, & Koratkar (1990) and Pogge & Peterson (1992) investigated the effect 
of variabihty on this relationship in C IV and Lya. The results can be summarized as follows: 

1. There is a global or object-to-object relationship with slope a 0.83 for CiV and 
a ^ 0.88 for Lya. 

2. There is a rather flatter intrinsic relationship with a ~ 0.4 that describes how the line 
and continuum relationship changes due to variability. At least some of the scatter in 
the global relationship is attributable to variability and the flatter slope of the intrinsic 
relationship. Moreover, the scatter in this intrinsic relationship is largely attributable 
to light-travel time effects within the BLR (i.e., changes in the line flux follow those in 
the continuum flux with some measureable time delay) . 

While the Baldwin Effect is well-established in most of the strong UV lines, any global 
Baldwin Effect in the Balmer lines seems to be quite weak (e.g., Dietrich et al. 2002). 
However, to the best of our knowledge, the question of whether or not there is an intrinsic 
relationship in the Balmer lines has not been speciflcally addressed. Our goal in this paper 
is to remedy this. We adopt an approach similar to that of Pogge & Peterson (1992) in 
attempting to remove light travel-time effects when determining the slope of the relationship. 
We will make use of numerous well-sampled optical spectra of NGC 5548 obtained by the 
International AGN Watch in this analysis to show that an intrinsic Baldwin Effect is present 
in Rp. 

2. Data Analysis 
2.1. The Data 

For this investigation, we used ten years of spectra from the International AGN Watch 
collaboration (Peterson et al. 2002 and references therein). These data are fairly inhomoge- 
neous, with differences in spectral resolution, quality, and amounts of contaminating light 
from the host galaxy. In order to minimize the effects of these problems, we restricted our 
analysis to the high-quality homogeneous data obtained with the Ohio State CCD spec- 
trograph on the Perkins 1.8-m telescope at the Lowell Observatory in Flagstaff, Arizona, 
between 1988 and 1998. The optical continuum and H/3 emission-hne fluxes for this subset 
are shown in Fig. 1. The continuum flux was measured in a 10 A wide bin at a wavelength 
of 5100 A in the rest frame of the galaxy, which has a redshift of 2; = 0.017. The continuum 
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underneath the H/3 line was modeled by linear interpolation between two windows on either 
side of the H/9 line, at wavelengths of 4710 A and 5100 A in the rest frame of the galaxy. 
The H/9 flux was then measured by integrating over the flux above this continuum. The 
measurement uncertainties in these spectral measurements were estimated by examining the 
differences in fluxes measured on very short time scales (i.e. less than or equal to 2 days). We 
assumed that flux differences on such short timescales are stochastic errors rather than real 
differences. We estimated the fractional error in the integrated continuum and hne fluxes to 
be 2.5%, consistent with previous estimates of the error in these types of calculations (e.g., 
Peterson et al. 2002). 



2.2. Light Travel-Time EfTects 

In this investigation, it is important to take into account the effects of light-travel time. 
Continuum radiation emitted by the central continuum source (presumably an accretion disk) 
is absorbed and reprocessed into emission-line radiation by diffuse BLR gas. Due to light 
travel-time effects within the BLR, we observe a time delay of order days to weeks between 
variations in the continuum and the response of the emission lines. A simple linear model 
(Blandford & McKee 1982) to describe the emission-hne response to continuum variations 
can be written as 

L{t)^ J ^{T)C{t-T)dr, (2) 

where C{t) and L{t) are the continuum and emission-line light curves, respectively, and r is 
the time delay, and ^(r) is the transfer function. The transfer function depends on the size, 
geometry, and structure of the BLR. Under typical conditions, the centroid of the transfer 
function is essentially the mean response time of an emission line to changes in the continuum. 
The mean response time, or lag, for an emission line can be found by cross-correlation of the 
light curves. 

In the case of NGC 5548, multiple years of monitoring optical continuum and emission- 
line variability has shown that the lljS time lag changes over time. More precisely, the lag 
varies with the mean continuum level in a fashion that clearly suggests that the response 
is dominated by gas in which physical conditions arc optimal for the H/? line. Thus, the 
mean response time for a given line is a dynamic quantity and eq. (2) is an overly simplistic 
model of the line response. However, it seems to be sufficient to account for time-dependent 
changes in the lag by taking yearly averages (Peterson et al. 2002). 

We characterize the time delay by the value of the centroid of the optical continuum/H/? 
emission-line cross-correlation functions for each year. The cross-correlation centroids arc 
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taken from Table 8 of Peterson et al. (2002); values range from a low of 11.37 days in 1992 
to a high of 26.44 days in 1998. For each point in the H/3 light curve, a corresponding lag- 
corrected continuum point was determined. The lag was subtracted from the time of each 
H/3 measurement to identify the time at which the optical continuum was best-correlated 
with the H/3 emission-line flux. The continuum flux at this time was determined by linear 
interpolation between the actual observations immediately preceding and following this time. 
Errors on interpolated continuum points were estimated by adding the errors of the preceding 
and following continuum points in quadrature, which probably slightly overestimates the 
actual uncertainty. 

2.3. Fitting Procedure and Results 

In Fig. 2, we show the total H/3 flux as a function of optical continuum flux, corrected 
for time delay as described in the previous section. The relationship is clearly non-linear. 
However, we have not yet accounted for constant spectral contaminants that have to be 
removed from the observed optical continuum and emission-line fluxes. The two major 
contaminants to the spectrum of importance to our investigation are (a) the contribution of 
the host galaxy to the optical continuum and (b) the H/3 narrow-line contribution to the total 
H/? flux. Assuming that these are the only two necessary corrections, we can write eq. (1) in 
terms of the obscrvables, namely the total continuum flux Fa(5100A) (i.e., AGN plus host 
galaxy) and the total B.f3 flux F(H/5totai) (i-e., broad component plus narrow component) as 

i^(H/?broad) = i^(HAotal) - i^(H/?narrow) = A(Fa(5100 A) 

-Pgalaxy) ■ (3) 

We can approach eq. (3) in two ways, flrst by a complete non-linear least squares treat- 
ment, and second by flxing the host galaxy and narrow- component fluxes, which are well- 
constrained by independent analyses, which reduces the problem to a linear least-squares 
exercise. We consider both of these in turn. 

Non-Linear Least Squares Fit. We employed a Marquardt algorithm as described by Press 
et al. (1989) to determine the constant components F(H/3narrow) and Fgaiaxy, the scaling factor 
A, and the Baldwin Effect index a, all as defined in eq. (3). The values of these parameters 
are shown in Table 1, along with uncertainties estimated by bootstrap resampling. The 
best-fit locus is shown in Fig. 2. 

Simple Power-Law Fit. We made use of the fact that two of the parameters in eq. (3) can 
be accurately determined by independent means. The host galaxy fiux at 5100 A is estimated 
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to be -Fgaiaxy = (3.37 ± 0.54) x 10~^^ ergs cm~^ A"^ (Romanishin et al. 1995) by using a 
model of the host-galaxy surface brightness distribution obtained by careful deconvolution 
of the nuclear point source and the host galaxy. 

We used a value for the narrow H/5 hne component of -F(H/3narrow) = 6.7 x lO'^^ergs 
cm~^ which was estimated by trial-and-error scaling of the [O ill] A5007 emission-line 

profile and subtracting the scaled profile at the position of the H/3 narrow line. This works 
best on high-quality data obtained during a low-fiux state when the narrow component of H/3 
is relatively prominent. A plausible narrow-H/3 model was found by scaling the [Olll] A5007 
emission line by 0.12 (±0.01), as shown in Fig. 3. After subtraction of these constant 
components, eq. (3) reduces to a simple linear least squares problem 

logF(H/3broad) = log A + alog [Fa(5100 A) 

-^galaxy ]■ (4) 

The linear regression was performed using the BCES (bivariate correlated errors and intrin- 
sic scatter) program described by Akritas & Bershady (1996). This technique accounts for 
measurement errors in both parameters. The results for BCES(Y||X) regression are given in 
Table 1, again with uncertainties in the fitted parameters determined by bootstrap resam- 
phng. The best fit function is shown in Fig. 4. 

3. Discussion 

Both methods of fitting show a clear Baldwin Effect, with a ^ 0.65 for the simple 
linearized fit and a ~ 0.53 for the non-linear least-squares fit, in both cases rather steeper 
than the values found by Pogge & Peterson (1992) for CiV and Lyct. The fit shown in 
Fig. 4 raises concern about either the parameterization of the relationship or the estimate 
of the background components (narrow-line H/3 and host-galaxy); in other words, if the 
independent estimates of the background components are correct, then as shown in Fig. 
4, most of the low-fiux points fall below the best-fit line, indicating that the slope of the 
relationship must be different in this regime. While the non-linear least-squares fit (Fig. 2) 
does not have this problem at the low-flux end, the background estimates from the fit are 
in both cases much too high, although certainly the constraint on the host-galaxy fiux is 
less stringent than that on the H/3 narrow component. To pursue this, if we adjust upwards 
from the Romanishin et al. (1995) estimate of the host-galaxy fiux, we find the best overall 
fit (as measured by x^) is achieved for a value of 4.6 x 10~^^ ergs s~^ cm' -2A-i(yield mg 
a ^ 0.52), which is, not surprisingly, close to the host-galaxy estimate we obtained (Table 
1) by the non-linear least-squares fit to eq. (3), the only difference being whether or not the 
narrow- line contribution is a free parameter or held fixed. A galaxy contribution of 4.6 x 
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10~^^ ergs s"^ cm^^ exceeds the Romanishin et al. (1995) estimate by ~ 2(7, but it is still 
(although only slightly) below the minimum published optical continuum flux (Peterson et al. 
2002). However, during the 2002 observing season, NGC 5548 went into the lowest-flux state 
recorded over the last 30 years, reaching a minimum 5100 A flux of (3.94 ±0.10) x 10~^^ergs 
s~^cm~^A~^ on JD2452471 (Sergeev 2002, private communication). Given this result, we 
doubt that the estimate of Romanishin et al. (1995) signiflcantly underestimates the host- 
galaxy contribution; it seems more likely that our parameterization is oversimplified. 

While Fig. 4 suggests that the relationship between the continuum and the H/? emission 
line changes at low flux levels, it is also possible that the relationship simply varies somewhat 
with time. As all of the low-flux points below the best-flt lines are from a single epoch (the 
fourth year of this monitoring program, during the period JD 2448713 - 2448848), we are 
unable to discriminate between flux dependence and time dependence. A complete analysis 
of the most recent faint-state data should clarify this. 

Finally, we as show in Fig. 3, our best estimate of the H/3 narrow component, based on 
a profile decomposition, is 6.7 x 10~^^ergs s~^ cm~^. This is much lower than the non- linear 
least-squares result of 1.8 x 10~^^ergs s~^cm~^ (Table 1; also see Fig. 2). It is very unlikely 
that the H/3 narrow component based on profile decomposition is underestimated. However, 
the non-linear least-squares result may include other less- variable contaminants to the broad 
H/3 profile, and this is under investigation. 

Depending on how we account for the constant components, we find possible power-law 
indices (eq. 3) in the range 0.52 < a < 0.65. In any case, we see clear evidence for an intrinsic 
Baldwin Effect in H/3, i.e., the line flux does not increase linearly with the continuum flux. 
To explore the physical nature of this relationship, however, we need to relate the H/5 flux 
to the ionizing continuum flux rather than the optical continuum flux, since the former 
actually drives the emission-line variations. While the ionizing flux (A < 912 A) is not 
directly observable, measurements closer to the ionizing continuum, at 1350 A, are available 
from the International Ultraviolet Explorer archive. Assuming a relationship of the form 

i^opt oc F^^, (5) 

Peterson et al. (2002) find (3 f» 0.56, based on contemporaneous pairs of 1350 A (Fuv) and 
5100 A (-Fopt) continuum fiuxes. However, by using the BCES regression, we found a slightly 
steeper slope, (3 — 0.67, which we will use in this discussion. We also note that in the case 
of NGC 5548, Galactic extinction is a small enough effect that it makes no difference in the 
slope /3 of eq. (5). 

Our analysis thus far has been based only on optical data, and we have assumed specif- 
ically that F(H/3) oc F^^. Combining this with eq. (5). we can write the relationship 
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between the H/3 emission line flux and the UV ionizing continuum as F(H/3) oc Fyv' ^^^^ 
0.35 <aP< 0.44. 

4. The Effect of Geometric Dilution 

The emission-line light curve is a convolution of the continuum light curve with a transfer 
function, as shown in eq. (2). The effect of the transfer function is to shift and broaden the 
emission-line light curve in time relative to the continuum curve; the effect of making the 
BLR larger is not only to shift the cross-correlation lag to larger values, but also to smear 
out the response in time, reducing the apparent amplitude of the response at any given 
time. We therefore consider here how important an effect this geometrical dilution of the 
line response might be in this particular case. We attempt to estimate this by convolving the 
real optical continuum light curve with a transfer function for a thin spherical shell, since 
this is a rectangular function that will have a relatively large dilution of the line response 
through time-smearing effects, and repeating the analysis described above. 

The effect of geometric dilution on the amplitude of variability of emission lines was 
explored by convolving the NGC 5548 continuum with thin-shell transfer functions of radius 
between 1 and 10 light days. We then used the original continuum data and the simulated 
H/3 emission to calculate the slope of the Baldwin Effect as discussed above. We found that 
the slope does decrease slightly with increasing geometric dilution; however, it does so at a 
very small rate. The difference in a is found to be only about 0.05 for rectangular transfer 
functions of half width 1 and 10 days. We tested transfer functions for other geometries, 
and found that, as expected, the effects of smearing were even smaller for narrower profile 
transfer functions. 



5. Comparison with Theoretical Predictions 

We can compare our empirical results on the intrinsic Baldwin Effect in B.f3 with the 
theoretical predictions made in a series of photoionization calculations by Korista et al. (1997; 
see their Fig. 3g). The quantity $(H) is defined as the number of ionizing photons arriving at 
the inner face of a BLR cloud per square cm per sec. If the shape of the ionizing continuum 
is fixed, the emission-line spectrum will be a function of $(H), which is proportional to 
the photoionization rate, and of the particle density n (cm~^), which is proportional to the 
recombination rate, as well as chemical abundances and cloud column density. Korista et 
al. show the strength of various emission lines normalized by the continuum flux (i.e., the 
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line equivalent widths) for various values of $(H) and n for different column densities and 
standard assumptions about elemental abundances. For NGC 5548, log$(H) = 20 at a 
distance of 12.6 light days from the central source (Korista & Goad 2001). The equivalent 
width used by Korista ct al. is defined as the ratio of the H/3 fiux to the continuum flux at 
1215 A. Using the above relationships, the equivalent width should be 

W^^c. FSi-\ (6) 

Using our values of a and we flnd that we expect the equivalent width to be pro- 
portional to the UV flux to the —0.56 or —0.65 power, depending how we account for 
the constant components. In other words, an order of magnitude change in $(H) should 
change the equivalent width by around —0.60 dex. Prom Fig. 3g of Korista et al. (1997), 
one can see that for log$(H) 20 and for hydrogen particle densities n(H)(cm~^) in the 
range 11 ^ logn(H) ^ 13, the slope of the log 1^/ log $(H) relationship is consistent with 
our results. The grids^ were investigated to find the slope of the relationship for varying 
values of column density and continuum shape with solar metallicity. Since the value of 
log$(H) pa 19.6 at a radius of 20 light days, where typically the H/5 response is maximized 
in NGC 5548, we determined the slope of the grid for varying values of the hydrogen particle 
density by linear interpolation along the grid. We found that regardless of the model, the 
slope of the logl^/log$(H) relationship reaches a maximum in the range —0.57 to —0.60 
when the value of the hydrogen particle density is in the range 12 ^ logn(H) ^ 12.5. 



6. Conclusions 

We have investigated the relationship between the variable H/3 emission line and optical 
continuum of NGC 5548. We discovered that an intrinsic Baldwin Effect does exist; after 
for correcting for hght travel-time effects, we flnd that the amplitude of broad H/3 variation 
is less than the amplitude of variation of the underlying optical continuum. By combining 
these results with the relationship between the amplitudes of variability of the optical and 
UV continuum, we determine the relationship between the HjS emission-line fiux and the 
UV continuum fiux, and find that our result is in good agreement with the predictions of 
photoionization calculations, even accounting for smearing of the line response by light travel- 
time effects within the H/3-emitting region. We suggest that similar analysis of multiple lines 
have potential as diagnostics of the physical nature of the BLR. 



^The photoionization grids calculated by Korista et al. (1997) can be found on the web at 
http://www.pa.uky.edu/'~korista/grids/grids.html. 
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Fig. 1. — The optical continuum (upper panel) and H/? (lower panel) light curves of NGC 
5548 from 1989 December to 1998 June. The continuum fluxes are in units of 10~^^ergs 
s^^cm""^ A~^, and the hne fluxes are in units of 10~^^ergs s~^cm~^. The dotted horizontal 
lines show the independent estimates of the constant components, specifically the continuum 
contribution from the host galaxy starlight (Romanishin et al. 1995) in the upper panel, and 
the narrow-line H/3 contribution in the lower panel. Flux measurements are in the observer's 
reference frame and are uncorrected for Galactic extinction. 
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Fig. 2. — H/3 vs. optical continuum flux for NGC 5548. The solid line shows the best fit 
to these data using a simple non-hnear form (eq. 3). The dashed lines show the constant 

components (host galaxy flux and H/9 narrow-line contribution) based on this non-linear fit. 
The dotted lines show the estimates of these quantities based on independent estimates, with 
the shaded area showing the la uncertainty in the host-galaxy flux. 
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Fig. 3.— The H/3 spectral region of NGC 5548 in a low state on JD2448726. The effect of 
attempting to remove the constant narrow-line component of H/3 is shown. The narrow-line 
H/? component was removed by using the [O ill] A5007 line as a template, scaled in flux by a 
factor of 0.12. 
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Fig. 4. — Broad H/3 flux as a function of optical non-stellar continuum flux for NGC 5548. 
The broad-line B.f3 flux is obtained by subtracting an H/9 narrow-line contribution of 6.7 x 

10^^"^ ergs s^^ cm^^ from the total line measurement. The non-stellar continuum flux is 
obtained by subtracting a starlight contribution of 3.37 x 10^^^ ergs s~^cm^^A^^ from the 
optical continuum measurements. The solid line is the best fit of eq. (4) to these data, with 
slope a = 0.648 (Table 2). 
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Table 1. Fits to the Optical Continuum-H/? Relationship 



Fit -f^(H/3narrow) -^galaxy Q; 

(1) (2) (3) (4) (5) 



Non-linear least squares 1.86 ± 0.28 4.85 ± 0.08 2.53 ± 0.17 0.526 ± 0.022 

(eq. 3) 

Power-law fit 0.67^ 3.37^= 2.18 ±0.09 0.648 ± 0.022 
(eq. 4) 



'^Units of 10~^^ ergs s~^ cm~^. 
•^Units of 10~^^ ergs s~^ cm~^ 
'^Parameter held fixed. 



